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ABSTRACT 

Hot stars are sources of X-ray emission originating in their winds. Although hydrodynamical simulations that are able 
to predict this X-ray emission are available, the inclusion of X-rays in stationary wind models is usually based on 
simplifying approximations. To improve this, we use results from time-dependent hydrodynamical simulations of the 
line-driven wind instability (seeded by the base perturbation) to derive the analytical approximation of X-ray emission in 
the stellar wind. We use this approximation in our non-LTE wind models and find that an improved inclusion of X-rays 
leads to a better agreement between model ionization fractions and those derived from observations. Furthermore, the 
slope of the Lx — L relation is in better agreement with observations, however the X-ray luminosity is underestimated 
by a factor of three. We propose a possible solution for this discrepancy. 
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1. Introduction 

Hot-star winds have been traditionally modelled assuming 
a spherically symmetric, stationary outflow. These assump- 
tions provide a convenient basis for studying different phe- 
nomena that may influence the wind structure. The line 
profiles and wind parameters predicted in this way are in 
good agreement with those found from observations (e.g., 
iPauldrach etalll2001l IVink et alJl200lLlKrtickall200^ . 

However, a stationary approach to hot-star wind mod- 
elling was for a long ti me known to be not fully adequate 
(|Lucv fe Solomonlll970D . Hydrodynamical simulations show 
the growth of strong shocks in the wind due to the so- 
called line-driven instability inherent to ra diative driving 
(|Owocki et al.lll988l IMdmeier et al.lll997al) . On the other 
hand, from the observational point of view the inhomo- 
geneities do not imprint a clear signature in hot-star wind 
spectra in ultraviolet and visible spectral regions, and only 
a detailed sp ectral analysis revea ls its possible pr esence 
(|Bouret et al.ll2003l IMartins et all [20051 iPuls et al.ll2006l) . 
Consequently, the influence of inhomogeneities (or clump- 
ing) on wind spectra has long been neglected. At present 
it is not clear whether the neglect of wind clumps in the 
formation of wind line profiles and infrared continua leads 
to an over estimate of mass-l oss rates derived from observa- 
tions (see IPuls et all l2008al for a discussion of this prob- 
lem) . 

The only clear observable signature of intrinsic wind 
non-stationarity is the existence of X-ray emission. This 
emission is directly observable for nearb y stars (e.g., 
iBerghofer et~aT] 119971 lAntokhin et~aT1 l2008f) and its exis- 
tence can be inferred from the influence on the ultraviolet 
spectrum for those stars for which a direct observation of 



their X-ray emission is not available (e.g., IBianchi et al.l 
l200l . 



The line-driven instability is not the only possible pro- 
duction mechanism of X-rays in hot star winds. Because the 
winds are highly supersonic, any mechanism which causes 
wind stream collisions may also lead to X-ray production. 
Consequently, in hot star binaries the X-rays may origi- 
nate due to the collisions of winds from individ u al stars 
fe-g.. iPrilutskii fe Usovl ll97oi iCooke et alJ fl978l iPittardl 
2009). The collisions of individual wind streams chan- 
neled by the magnetic field may also cause X-ray em ission 
(jBabel fc Montmerlel[l997l lud-Doula fe Owockill2002D . The 
former mechanism is important in binaries and the latter 
one in stars with a sufficiently strong magnetic field. In 
this paper we concentrate on a more general mechanism, 
operating in all radiatively driven stellar winds, i.e., the 
line-driven instability. 



A central observational finding related to the X-ray 
emission is the dependence of the X-ray luminosity Lx 
on the total lumi nosity L via the appro ximate relation 
L x ~ 10~ 7 L (e.g.. IChlebowski et al.lll989D . This relation- 
ship is still not explained by wind theory. A s sumin g a con- 
stant X-ray filling factor, lOwocki fe Cohen] (|1999( ) showed 
that for stars with optically thick winds the observed X- 
ray luminosity scales with the mass-loss rate as Lx ~ M, 
and for stars with optically thin winds as Lx ~ M 2 . As 
the predicted wind mass-loss rate scales with the stellar 
luminosity as M ~ £ 1/q ' (|Kudritzki fc Pulsll2000l ). where 
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a' k 0.6 for luminuou s O star^]] (e.g., Vink et all l2000t 
IKrticka fe Kubatl 120091 iPuls et all f2008bf) . the predicted 
slope of the Lx — L relation is steeper than the observed 
one. This difference may origina te in the radial depen- 
dence of the X-ray filling factor (|Owocki fc Cohen| [l999) 
or may be caused by t he dependence of the co oling length 
on the wind density (Kr ticka fc Ku bat 200 l|). Note also 
that macroclumping (|Feldmeier et al.ll2003l lOskinova et all 
12004 lOwocki et al.ll2004 lOskinova et al.ll2006D may lead to 
a decrease in the effective opacity in the X-ray region and 
affect the X-ray luminosity. 

Although a real hot-star wind should be far from sta- 
tionary and spherically symmetric (on small scales), the 
success of wind modelling applying these approximations 
motivates us to incorporate the inhomogeneities found in 
time-dependent hydrodynamical simulations into a spher- 
ically symmetric, stationary wind model. These models 
can serve as an efficient tool to study hot-star winds un- 
til more elaborate hydrodynamical simulations that consis- 
tently take into account the non-LTE effects become avail- 
able. 

There have been earlier attempts to include X-ray 
emission in non-LTE wind models. T hey were either 
based on simplified analy tical models (|Feldmeier et al.l 
Il997bl iMartins et all [200l . or the X -ray emission was 
included usi ng ad-hoc free pa ramet ers (aka the "filling 
factor", e.g.^ MacFarlaneeLalJ 11994 IMartins et all 120051 
IKrticka fc Ku bat 2009) describing the hot wind part. Here 
w e use the results o f availa ble hydrodynamical simulations 
of iFeldmeier et aTl (|1997af) to directly describe the X-ray 
emission in a compact form and include it in our non-LTE 
wind models. 



2. Hydrodynamical simulations 

The hydrodynamical simulations we employ here to esti- 
mate the X-ray emission from O stars were performed us- 
ing the smooth source function method (SSF, Owocki 1991, 
lOwocki fc PulsllT996l ). This corresponds roughly to a for- 
mal integral of the radiation force using a pre-specified 
line source function, here the optically thin source func- 
tion, S ~ r~ 2 . The line-driven instability that is responsi- 
ble for the formation of X-ray emitting shocks in the wind 
is captured by a careful integration of the line-optical depth 
with a resolution of three observer-frame frequency points 
per line Doppler width and accounting for non-local cou- 
plings in the non-m onotonic velocity law of the wind (see 
lOwocki et al.l ll988. for details). The angle integration in the 
radiative flux is limited to one single ray which, however, 
is not the radial ray but hits the star at « 0.7i?». This ac- 
counts with sufficient accurac y for the finite-disk correction 
factor (|Pauldrach et al.lll986l ) and avoids the critical-poin t 
degeneracy of the point-star CAK model (|Poe et al.H l990V 
To avoid any artificial overestimate of the X-ray pro- 
duction in the wind, special care was taken in the nu- 
merical calculation of the line force on the staggered spa- 
tial mesh in ord er to correctly reproduce the line-drag ef- 
fect ()Lucvl[l984K which is responsible for a partial reduc- 
tion of the instability, esp ecially in layers close to the star 
(jOwocki fc Rvbickil Il985h . The pure hydrodynamics part 



1 Note that the situation may be different in late O and B 
supergiants, for which the value a' ~ 0.8 — 1 might be a more 
appropriate one (Bcnaglia et al. 2 0071 ). 



of the code is a standard van Leer solver following stan- 
dard prescriptions for such a scheme: staggered mesh; op- 
erator splitting of advection and source ter ms; ad vection 
terms in a conservative form using van Leer 's lj!977l ) mono- 
tonic derivative as an optimised compromise between sta- 
bility and accuracy; Richtmyer artificial viscosity; and non- 
reflecting b oundary conditions (|Hedstromlll979LlThompson1 
Il987l . Il990l ). The numerical collapse of post-shock radia- 
tive cooling zones caused by th e strong oscillatory thermal 
instability (jLanger et al.lfl98l[ ) is prevented by artificially 
modifying the radiative cooling function below a certain 
temperature at which X-ray emissivity is still small. Further 
details can be found in iFeldmeierl (|l995l ). 

As seed perturbations for unstable growth we introduce 
a turbulent variation of the velocity at the wind base, at a 
level of roughly one third of the sound speed. This perturba- 
tion is obtain ed from a sim ple quadrature of the Langevin 
equation (see Risken 1996) with a coherence time in the 
friction term of 5000 seconds, which is well below but not 
too far off the acoustic cutoff period of the star at which, 
according to the theorem by Poincare (see iLambl Il945l ). 
acoustic perturbations of the photosphere should accumu- 
late. The power spectrum E(k) of this Langevin turbulence 
has a spectral index of —2, which is close to the Kolmogorov 
index of —5/3 for the universal, inertia subrange of turbu- 
lence. 

The variety of dynamical structures in the instability- 
induced, line-driven wind turbulence (which grows out of 
but is not identical to the turbulence applied in the photo- 
sphere via boundary conditions) is still largely unexplored, 
but seems to be similarly intricate to that found in other 
turbulent flows in ordinary fluids and gases or in astro- 
physical MHD settings. There are indications of the pres- 
ence of a quasi-continuous hierarchy of density and velocity 
structures in the wind, similar to that foun d in supersoni c 
Burgers turbulence with shock cannibalism (|Burgerslll974 ). 
and there are also indications of a separation of dense struc- 
tures into two distinct families, which we address under the 
names "shells" and "clouds." The shells are formed close to 
the star, in a first stage of unstable growth. They do not 
collect all the wind material but rather one half of it, since 
the negative-velocity pert urbations remain unaff ected by 
the line-driven instability (|MacGregor et al"1ll979f h Further 
out in the wind, in a second stage of unstable growth, clouds 
are "ablated" from the outer rim of the remaining mass 
reservoir at CAK densities, and are accelerated by the stel- 
lar radiation field through the emptied regions, eventually 
colliding with the next-outer shell , and producing observ - 
able X-ray flashes in this collision (jFeldmcie r et al.lll997al ). 

Without photospheric Langevin turbulence, clouds do 
not occur and the mass reservoir is continuously fed into 
the next outer shell via a thin s tream of fast gas that by 
far cannot (see lHillier et al.ll 19931 ) account for the observed 
X-ray emission from O stars. To test for the relevance of 
the specific form of photospheric turbulence, we also calcu- 
lated models with a "stochastic" photospheric sound wave 
as seed perturbation for the instability, i.e. a sound wave 
with stochastic variations in period, amplitude, and coher- 
ence time, which resulted in largely the same results. 



3. X-ray emission from hydrodynamical simulations 

To incorporate the results of hydrodynamical simulations of 
£ Ori A wind in non-LTE wind code in a manageable way, 
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we approximate the emission from hydrodynamical simu- 
lations as a polynomial function. This could be done in 
two ways. The first way is to calculate the X-ray emission 
based on a local values of hydrodynamical variables, and 
then determine a function that best approximates its ra- 
dius and frequency dependence (see Sect. 13. ip . The second, 
somewhat simpler way, is to find a polynomial that fits the 
temperature structure of the simulation, and then use this 
polynomial to calculate the emission (see Sect. 13. 2| ). 

3.1. Spectrum fitting 

First, the wind structure from hydrodynamical simulations 
is used to calculate the X-ray emissivity as a function of 
radius and frequency. We divide the wind model into 29 
equally spaced radial bins and calculate the total wind X- 
ray emission in each radial bin at the timestep k as an 
integral 



l J k 



r 2 n c (r)n^(r)A(v.j,T(r)) dr, 



(1) 



where the radius denotes the radial boundary of the i- 
th bin, n c (r), nn(r) are the number densities of free elec- 
trons and hydrogen, and the X-ray emissivity A(i/j, T(r)) is 
calculated using the Raym ond-Smith X-ray spectral code 
(Raymond & Smith Il977l Raymond Il988l ) for uniformly 
spaced frequencies Vj with the temperature T(r). The vari- 
ables n c (r), nn(r ), and T(r) were taken f rom the simulation 
of C Ori A wind (jFeldmeier et al.lll997ah . We calculate the 
quantities 

1 I 



1L 

Ai/ rrii 



P.. — 



where 



r 2 n e (r)hn(r) dr, 



(2) 



(3) 



lij is the time average of and Av is the frequency spac- 
ing. n e (r) and hn(r) are the electron and hydrogen num- 
ber densities from a stationary wind model corresponding 
to hydrodynamical simulations. This stationary model is 
calculated by adopting the mean mass-loss rate from the 
hydrodynamical simulations (jFeldmeier et al.l [l997a) and 
using a standard beta velocity law. 

To provide a simpler description of the X-ray emission, 
we fit the values of £ij/10 _23 ergs -1 cm -3 keV" 1 by (see 
Fig. [J) 

log£(x, 7) = a(x) + b(x)-y + c(x)-y 2 + d(x)~f 3 , (4) 

where 7 = log(_E/lkeV) is the logarithm of energy E ex- 
pressed in keV, and x — r/R*. The radial variation of the 
polynomial coefficients is given by the expressions 

a(x) = -6.087 + 1.836x6 - 0.211^ - 0.112 (x - x 6 ) , 
b(x) = -4.124 + 0.232x 6 - 0.132 (x - x 6 ) , 
c(x) = 0.651 - 0.456a; 6 , 
d(x) = 2.457- 0.369a; 6 , 



(5) 



X6 



X, 

6. 



x < 6, 
x > 6. 



Note that £ — for x < 2. The X-ray emissivity in the 
non-LTE model can be then approximated by 



47T 



n e n H ^(?"/i?*,7) 10 



23 ergs ^m 3 keV 



where n e , riH are the number densities of free electrons and 
hydrogen in the non-LTE model. /* is the scaling factor 
(see below, Eq. (|12p). which allows one to use Eq. ([6]) also 
for stars other than £ Ori A. In Eq. we introduced the 
factor 1/47T because the emissivity A(v, T) is assumed to be 
that in all spatial directions. 

3.2. Temperature distribution fitting 

Instead of fitting the X-ray emissivity from simulations di- 
rectly, it is also possible to fit the temperature distribution 
of X-ray emitting gas. This approach does not explicitly de- 
pend on a particular functional form of the X-ray emissivity 
A(u, T), but is also not completely independent of it, as the 
distribution of temperatures in the numerical simulations 
depends on the adopted form of the cooling function. 

As the amount of emitted X-rays per unit volume of gas 
depends on the square of the density n 2 , it is not realistic 
to fit just the distribution of temperatures. To account for 
the correlation of the distribution of n 2 with temperature, 
we evaluate for each timestep k th e quantity similar to th e 
differential emission measure (e.g., Dzifcak ova et alj|2008h 



DEM 



ijk 



r 2 n c (r)nn{r)8j(T(r))dr, (7) 



where S 3 (T(r)) is equal to 1 for Tj - AT,-/2 < T{r) < 
Tj + ATj/2 and otherwise. Here we introduced the log- 
arithmically spaced grid of temperatures Tj with grid size 
ATj. The normalised DEM is introduced as 



1 DEMj 
X> ■ = - 

13 ATj m, 



(8) 



where DEM^- is the time average of DEMy^. 

To provide a simpler description of X>y we approximate 
it by a polynomial fit as 



log£>(a;, t) — a(x) + b(x) t, 



(9) 



(G) 



where t = log(T/l K) is the logarithm of temperature ex- 
pressed in units of Kelvin, and x = r/_R*. The radial vari- 
ation of the coefficients a(x), b(x) is given by 

a(x) = 3.383 + 1.509x6 - 0.279x 2 + 0.416 (x - x 6 ), 

b(x) = -2.435 + 0.0564x6 + 0.00724xg - 0.0783 (x - x 6 ). 

(10) 

Note that V = for x < 2. For given number densities of 
free electrons and hydrogen the X-ray emissivity is 

nx{r,v) = j-n c n H JJ V{r/R r ,\ogT)A„(T)dT, (11) 

where we set 7] = 10 5 K, T f = 2.5 x 10 7 K, and A U (T) is 
the X-ray emissivity per unit of frequency. A comparison 
of T> as derived from the simulations and our fit is given in 

Fig.m 

We checked that both approximations of X-ray emis- 
sion from hydrodynamical simulations (Eq. ((6|) or (|11[) ) and 
models based on exact data give similar results. 
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Fig. 1. Le/t: The frequency distribution of X-rays expressed in terms of log(^ij/10 23 ergs 1 cm 3 keV 1 ) as a function 
of radius and frequency (see Eq. @). Right: Fit to the X-ray emissivity expressed as log£ (Eq. (]3|)). 
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Fig. 2. Left: The temperature distribution of X-ray emitting gas expressed as log XV, , as a function of radius. Right: Fit 
to the temperature distribution expressed as log V. 



3.3. Scaling of X-ray emissivity with stellar parameters 

The X- ray emission in t he hyd rodynamical wind simula- 
tions of iFeldmeier et all (|1997a| ) was obtained for a mass- 
loss rate G ri A — 3 x 10~ 6 M Q /year, radius ori A = 
24 R , and terminal velocity ori A = 1850 kms -1 cor- 
responding to a hydrodynamical wind model of £ Ori A. 
Most of the X-rays used in these simulations originate due 
to the collisions of fast clouds with dense shell fragments, 
where the cloud formation is triggered by a turbulent pho- 
tospheric seed perturbation. Consequently it is not clear 
whether the derived analytical approximations are realis- 
tic also for other stars. To account for the difference be- 
tween the X-ray emiss ion in individual star w ind models 
and C Ori A model of IFeldmeier et all (|l997al ). we intro- 
duced a scaling factor /* in Eqs. ([6]), (fTTj) . 

The X-ray emissivity is given by the rate of energy dis- 
sipation by shocks. Neglecting possible dependencies on ve- 
locity and base perturbations, we expect the rate of energy 
dissipation to be proportional to the wind density. However, 
this is not what we would get from the scaling ?/x ~ p 2 in 
Eqs. ([6]), (TTT]) with neglected dependence on the individual 
wind parameters, i.e. with /* = 1. The reason is that we did 



not take into account that the fraction of X-ray emitting 
material may not be the same for all stars. In the winds with 
similar abundances, the radiative cooling is more efficient 
in stars with denser winds than in stars with low-density 
winds. In the winds with low density, a larger fraction of 
the wind gas may be involved in radiating away the shock 
thermal energy. 

Based on these arguments, we could introduce the scal- 
ing of /* with the density p at a given point in the wind 
as /* ~ 1/p. However, this is not the most convenient scal- 
ing, as it would introduce additional dependence of /* on 
the radius, which is already accounted for in Eqs. ([6]), (fTTj) . 
Consequently, for each star the value of /* should be fixed. 
Taking into account that p ~ M/(r 2 v), we introduce a 
representative wind density p — M / (R^Voo) and in our cal- 
culations use 



r 



PC Ori A _ M$ Ori A R* 



P 



M 



R C Ori A V < Or 



(12) 



where M, and are radius, wind mass-loss rate, and 
terminal velocity of individual stars. 
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Fig. 3. Left: The frequency distribution of X-rays derived using the "filling factor" approach expressed in terms of 
log(^j/10 -23 ergs" 1 cm" 3 keV" 1 ) as a function of radius and frequency (cf. with Fig. [J). Right: The temperature of 
X-ray emitting gas expressed as a function of radius (cf. with Fig. [21). 



4. Comparison with the "filling factor" approach 

It is informative to compare the X-ray spectrum and the 
temperature distribution function derived from hydrody- 
namic simulations in the previous section with those de- 
rived using simple "filling factor" approach. The filling fac- 
tor determines the amount of X-ray emitting material. 

We calculated the wind model of £ Ori A with X- 
ra y emission included afte r the "filling factor" approach 
of IKrticka fc KubaH J2009) and derived the resulting fre- 
quency distribution of emitted X-rays and the temperature 
distribution of X-ray emitting gas (see Fig. [3]) . The adopted 
model assumes that the temperature of the X-ray emitting 
gas is given by the Rankine-Hugoniot condition with the ve- 
locity difference related to the wind speed. Consequently 
the temperature of X-ray emitting gas in Fig. [3] increases 
with radius and X-rays are harder in outer wind parts. 

This is different from the results of hydrodynamical 
simulations. In the simulations the highest temperatures 
of X-ray emitting gas are achieved close to the star (see 
Fig. [5]). As a result, the most energetic X-rays are emit- 
ted close to the star (see Fig. [T|) . Such a temperature dis- 
tribution is supported by ob servations, as was inferred by 
Wal dron fc Cassinellil (|2007l ) from their extensive analysis 
of high-resolution X-ray spectra of O stars. Consequently, 
two important features of the X-ray emission emerge from 
our simulations (and from observations): the decrease of 
the temperature of X-ray emitting gas with radius in the 
outer wind, and the multitemperature distribution of X-ray 
emitting gas at a given point. 



5. Non-LTE models 

5.1. Model assumptions 

As an application of derived analytical formulae, we include 
X-ray emission parameterised by Eqs. (fTTj) and (jT2J) into our 
stationary, spherically s ymmetric non-LTE wind model (see 
IKrticka fc Kubatll2009l for a detailed description). The X- 
ray emission is included in the radiative transfer equation of 
our models. In the following we discuss just the temperature 
fitting approach as it seems to be more straightforward. The 



hydrogen and free electron densities in the X-ray source 
term are taken from the non-LTE models. 

The excitation and ionization state of the considered ele- 
ments is calculated from statistical equilibrium (non-LTE) 
equations. We use atomic data from the Opacity Project 
and the Iron Project (Seaton Il987l iFernley et'aTJ Il987l 
Luo & Pradhan [1989|, Sa wey fc Berrington [1992], Seaton 
ct al. 119921 Bu tler et al. 119931. Nahar & Pradhan 119931 
Hum mer et al. 19931 Bautista Il996l Na har fc Pradhan 
119961 Zh ang 119961 Bautista & Pradh an 1 19971 Zhang & 
Pradhan 119971 . Chen & Pradhan [19991 ) . A significant part 
of the ionic models was ta ken from the TLUSTY code 
(|Lanz fc Hubenvj[2003ll2007T h For phos phoru s we employed 
data described by Pauldrach et al. f2001). Auger pho- 
toionizatio n cross sections from indiv idual inn er-shells were 
taken from lVerner fc Yakovlevl(ll995L see also lVerner et al.l 
I1993T) - and Aug er yields were taken f rom iKaastra fc Mewe 
(|l993l ) . We use lAsplund et al.l (|2005l ) solar abundance de- 
terminations in the code. 

The radiative transfer equation is split in two parts, 
i.e., the continuum and line radiative transfer. The radia- 
tive transfer in the continuum is solved using the Fea utrier 
method i n sphe rical coordinates (Mihalas & Hummer |1974| 
or Kubat I1993D with inclusion of all free-free and bound- 
free transitions of the model ions. The radiative transfer in 
lines is solved in the Sobolev approximation (e.g., Castor 
Il974l ) neglecting continuum opacity and line overlaps. 

The radiative f orce is calculated in the Sobolev approxi- 
mation (see Castor IT9741) . The corresponding line data were 
extra cted in 2002 fr om th e VALD database (Piskunov et al. 
1995, Kupka et al. Il999f ). The radiative cooling and heat- 
ing terms are calculate d usin g the electron thermal balance 
method (Kubat et al., Il999f ). For the calculation of these 
terms we use occupation numbers derived from the statisti- 
cal equilibrium equations. Finally, the continuity equation, 
equation of motion, and energy equation are solved itera- 
tively to obtain the wind density, velocity, and temperature 
structure. 

The shortest wavelength considered in the models is 
4.1 A and the longest X-ray wavelength is defined as 100 A. 
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Table 1. Stellar parameters of studied O stars. 
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5.2. Test stars 

The purpose of our study is to investigate the trends in 
the X-ray emission known for all OB stars, rather than 
to study X-ray properties of individual objects. Therefore, 
we compiled a list of O stars that includes the objects of 
various luminosity classes, masses, and binarity status. By 
doing this, our list may be considered as representative of 
a diverse population of O stars found in real star clusters. 
This allows us to compare the trend in X-ray luminosity 
of our synthetic O star popula t ion w ith t he real clusters 
recent ly observed by I Sana et al.1 (|2006l ) and lAntokhin et all 
(|200l . 

T o address the "weak wind problem" (jBouret et al.1 
120031 iMartins et ail l2005h . we also included stars that 
should be subject to it, i.e., their observed wind lines are 
much weaker than expected (the bottom four stars below 
the horizontal line in Table [T]). 

are taken from 
Markova et al. 



The stella r para meters (see Table 
Lamer s et al. (119951) . Repolust e t al. (1200 



(|2004h . and lMartins et alj (120051). Note that the parameters 
deri ved by Repolu st _et L al. d200j), Markova et al. (|2004j) . 



and IMartins et alj (|2005l ) were obtained using blanketed 
model atmospheres, i.e., they should be more reliable than 
the older ones. Stellar masses were obtained using evolu- 
tionary tracks (Schaller et al. Il992[ ). The mass-loss rates 
in Table [T] were theoretically predicted using our models 
discussed in Sect. 15.11 
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Fig. 4. The dependence of the total X-ray luminosity on 
the bolometric luminosity calculated using non-LTE mod- 
els with X-ray emissivity after Eq. (fTTj) and Eq. ([12"]) (filled 
and empty circles) for individual stars. This is compared 
with the mean observational rel ations (corrected for th e in- 
terstellar absorpti on) derived bv lAntokhin et all (|2008T ) and 
ISana et ail (|2006l with uncertainty). Empty circles denote 
values for stars showing the "weak wind problem". 



6. Non-LTE models with X-ray emission from 
simulations 

The inclusion of X-ray emission does not significantly mod- 
ify the ionization fractions close to the star below the 
critical point radius where the mass-loss rate is fixed. 
Consequently, X-rays do not significantly modify the wind 
mass-loss rate, bu t may modify the wind terminal velocity 
by a few percent (|Krticka fc Kubatll2009l) . 



6.1. X-ray luminosity 

A comparison of the predicted X-ray luminosities for in- 
dividual stars and the mean observational trends is given 
in Fig. |4] The predicted X-ray luminosities for stars with 
optically thick winds (L > 1O 5 L ) are on average lower 
roughly by a factor of three than the obsery e d ones . This 
corresponds to the results of iFeldmeier et al] (|l997al) . 

The derived slope of the Lx — L relation for stars with 
optically thick winds, Lx ~ L , is in better agreement 
with observations. The cause of the decrease of the slope 
(compared to results obtained for a fixed filling factor as a 
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The reason for this behaviour is the inverse proportional- 
ity of X-ray opacity an d energy (e.g. lOskinova et al.l l2006. 
IKrticka fc Ku bat 2009). For stars with low luminosities, the 
wind density is low, and most of the emitted X-rays will be 
emergent. On the other hand, for stars with high luminosi- 
ties, density and opacity values become large, more flux is 
absorbed in the soft than in the medium energy band, con- 
sequently the emergent X-rays become harder. This also 
explains why the predicted hardness of X-rays is greater 
for giants and supergiants than for main sequence stars 
(see Fig. [5]). Finally, although from our models it follows 
that the predicted dispersion of X-ray luminosities is in- 
deed higher in the hard band than in the medium and soft 
energy bands, there is still a clearly discernible Lx — L re- 
lation even in the hard band. 



Fig. 5. Comparison between the predicted ratio of X-ray 
luminosities in the medium and soft energy bands for indi- 
vidual stars and the mean trend derived from observations. 
Main sequence stars are denoted with empty symbols, gi- 
ants and supergiants with filled symbols. 

free parameter, IKrticka fc~K ubat 200^) can be most easily 
understood from Fig. [I] which shows that £(x,j) decreases 
with x (i.e., with radius). Consequently, the filling factor 
also decreases wi t h rad ius, supporting the explanation of 
lOwocki fc Cohenl (|1999( ). Also the dependence of /* on the 
stellar parameters contributes to the decrease of the slope of 
the Lx — L relation, but the contribution of this parameter 
is sensitive to the stellar sample considered. 

For stars with optically thin wind (L < 10 5 L Q ), even 
the improved inclusion of X-ray emission is not capable 
of explaining the observed slope of the Lx — L relation. 
As indicated above, the dependence of the predicted X-ray 
luminosity of these stars on stellar luminosity is steeper 
than for stars with high luminosity. However, this is not 
supported by observations. The Lx — L relation of low lu- 
minosity stars is either the same as that of high luminosity 
stars , or, for B stars, the slope is even less steep (|Sana et al.l 
2006). The reason for this discrepancy between observation 
and theory is unclear. Note that for these stars the shock 
cooling le ngth could be comparable with the hydrody nam- 
ical scale (|Cohen et al.ll200a IKrticka fc Kubatl l2009f ). The 
fact that a substantial fraction of the wind is too hot to 
give a significant signature in the ultraviolet spectrum of 
the s tar could be an explan ation for the weak-wind prob- 
lem (|Krticka fc Kubatll2009D . 

6.2. The energy distribution of X-rays 

ISana et all (|2006fl and I Antokhin et all (|2008f) divide the X- 
rays into three energy bands, soft (0.5 — 1.0 keV), medium 
(1.0-2.5 keV), and hard (2.5 - lO.OkeV). The observations 
show that the slope of the Lx — L relation is roughly the 
same in the soft and medium band as in the total X-ray 
emission, whereas the observations in the hard band show 
a large dispersion. 

Our calculations show a somewhat different result, as 
can be seen from Fig. [5) where we plot the ratio of the 
X-ray luminosities in the medium (Lm) and soft (Lg) en- 
ergy bands. The predicted slope of the Lx — L relation in 
the medium band is steeper than that in the soft band, 
consequently the ratio of Lyi/Lg increases with luminosity. 



6.3. Ionization fractions 

The presence of X-rays may also influence the ionization 
fractions of highly ionised ions. We plot in Fig. [5] the pre- 
dicted ionization fractions in comparison with those derived 
from observation as a function of effective stellar tempera- 
ture. Note that here we use the io nization fractions derived 
for LMC stars (|Massa et al.ll2003[ ). because their observa- 
tional sample includes a large number of ionization stages, 
the fractions of which are given as a function of velocity, 
and the wind parameters of LMC stars are not significantly 
different from Galactic ones. Since several parameters ex- 
cept the effective temperature (e.g., the wind density) may 
influence the ionization fractions, the graphs in Fig. [6] are 
not monotonic. 

In Fig. [5] we also plot the results of IKrticka fc Kubatl 
(|2009fl . which were obtained using a simpler inclusion of 
X-rays assuming a constant filling factor (see also Sect.0]). 
Still, the emergent X-ray luminosities from these models 
roughly correspond to the observed ones. Rather surpris- 
ingly, although the non-LTE models calculated with X-ray 
emission from hydrodynamical wind simulations give a too 
low emergent X-ray luminosity, the ionization structure of 
these models corresponds in general much better to the 
trends derived from observations. 

This better agreement between theoretical and obser- 
vational ionization fractions may be due to an overall 
lower X-ray emissivity. To test this, we calculated mod- 
els with a simpl i fied t reatment of X-ray emission after 
IKrticka fc Kubatl ([20091 ) . but assuming the same X-ray lu- 
minosity as that predicted by hydrodynamical simulations. 
Our results show that even in this case the use of X-ray 
emissivity based on hydrodynamical simulation gives bet- 
ter agreement with observations. We conclude that the hy- 
drodynamical simulations give a more realistic dependence 
of the X-ray emissivity on frequency. 

7. Discussion 

The remaining discrepancy between theory and observa- 
tion may originate in our neglect of macroclumping (some- 
times also termed porosity) . The macroclumping may cause 
a lower opacity in the X-ray domain, which leads to 
a higher predi c ted X -ray luminosity. From Fig. (17) in 
lOskinova et al.l {2004), it follows that an increase of X- 
ray luminosity by a factor of 3 would be easily achievable. 
Accounting for macroclumping reduces the wavelength de- 
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Fig. 6. Ionization fractions as a function of the effective temperature for individual stars from our sample (only stars 
with M > 10 M Q /year are included here) at the point where the radial velocity v — 0.5uoo. Filled circles • refer to the 
present models with X-ray emission from hydrodynamical simulations , open circles denote values taken from former 
non-LTE models with X-ray emission descri bed using fill i ng fa ctor (jKrticka fc Kubatll2009fh The ionization fractions 
derived from observations were adopted from lMassa et alj (|2003l for LMC stars, plus signs +). 



pendence of opacity. In its limit ( when clumps ar e fully 
opaque), the opacity becomes grey (lOskinova et alj |2004). 

Wind inhomogeneities should also affect the ioniza- 
tion fractions. Optically thin inhomogeneities (often re- 
ferred to as "m icroclumping" ) lead generally to lower ion - 
ization stages (|de Koter et all 120081 jKrticka et all 120081) . 
hence clumping with increased X-ray emissivity could be 
another way to explain the X-ray luminosities and ioniza- 
tion fractions derived from observation. 

Our hydrodynamical wind simulations assume 1-D 
spherical symmetry, thus all wind structures correspond to 
full spherical shells. Therefore, we have at present a very 
limited knowledge of real wind clumping. 

To be consistent with the hydrodynamical simulations, 
we use th e Raym ond-Smith X-ray spectral code (Raymond 
& Smith 119771 Raymond |1988[ ). although it would be 
possible to use codes based on up-to-date atomic data. 
Our test using the APEC subrou tine available in XSPEC 
(jArnaudl 119961 ISmith et all 120011 ) showed that there is a 
relatively good agreement between the fits using this and 
the Raymond-Smith code Eq. ([5]). 

8. Summary 

We provide analytical approximations of the X-ray emis- 
sion pred icted by hydrodynamica l simulations of hot star 
winds bv lFeldmeier et all (|1997aD . The X -ray emission de- 
rived from the hydrodynamical simulations has two distinc- 
tive aspects. First, the temperature of X-ray emitting gas 
decreases with radius in the outer wind. Second, the tem- 
perature of the X-ray emitting gas is described by a distri- 
bution function, which is more realistic than assuming just 
one temperature at a given point. 

We include these ap proximations in non-LTE models 
(jKrticka fc Ku bat 2009) and for selected stars compare the 
resulting X-ray luminosities, energy distribution of emer- 
gent X-rays, and ionization structure with observational 
results. We conclude that although the predicted X-ray lu- 
minosity is slightly underestimated, the resulting ionization 
structure is in relatively good agreement with observations. 
Earlier papers debated whether the theory of line-driven 



winds can explain the observed Lx ~ L relation. Our re- 
sults reproduce this scaling. 
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